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ABSTRACT In this paper, an attribute control chart using the repetitive sampling under the neutrosophic 
statistics system is discussed. The necessary measures of the proposed control chart under the neutrosophic 
statistics system are given. The control chart coefficients of the proposed control chart are determined using 
an algorithm under neutrosophic statistics system. The efficiency of the proposed control chart in terms 
of neutrosophic average run length is discussed over the existing control chart under neutrosophic statistics 
system. From the comparison studies, it is found that the proposed control chart under neutrosophic statistics 
system is more sensitive in detecting a shift in the process as compared with the existing control chart under 
neutrosophic statistics system. An industrial application of the proposed control chart under neutrosophic 


statistics system is also given. 


INDEX TERMS Attribute chart, average run length, control chart, neutrosophic logic, shift. 


I. INTRODUCTION 

The Shewhart control charts have been widely used in the 
industry to monitor the controlling the number of non- 
conforming items or defective items during the manufactur- 
ing process. The Shewhart control chart consists of three 
control limits are known as upper control limit (UCL), lower 
control limit (LCL) and the central limit (CL). A process is 
declared as out-of-control if the plotting statistic exceed the 
LCL or UCL. The control charts help to keep the process 
near the CL and minimize the non-conforming items. The 
Shewhart X-bar control chart is applied when the data is 
continuous and Shewhart np control chart is used when data 
is discrete. The X-bar control chart is more informative than 
the attribute control chart. However, the earlier control chart 
cannot be applied when the purpose is monitoring of non- 
conforming items. In addition, the np-control chart is quite 
simple and easy to apply in the industry. Further, as men- 
tioned by [1], the attribute control chart can also be applied 
when the quality of interest is categorical. However, fuzzy 
approach [2] has the ability to monitor more than one quality 
characteristics simultaneously. 

The fuzzy logic is applied when the parameters or the 
observations are imprecise or unclear. The fuzzy attribute 
control charts can be applied to monitor the categorical char- 
acteristics in indeterminacy environment. According to [3] 
“Due to the emphasis on the user's feelings and psychological 


factors, there are many fuzzy attributes of quality from the 
fitness-for-use viewpoint. As a result, there are not only two 
distinct judgments (applicable or inapplicable) when evalu- 
ating the quality from a fitness point of view. In this sense, 
considering the fuzzy property of fitness-for-use quality is 
more practical". The fuzzy chart for monitoring attribute 
data was proposed by [4]. The economic aspects of fuzzy 
attribute chart were considered by [5]. The application of 
such fuzzy charts is given by [6] and [7]. Later on, [8], [9], 
and [10] designed charts to monitor imprecise data. The 
fuzzy charts using weighted average methods was considered 
by [2] and [11]. The range chart was proposed by [12]. 
A detailed discussion and applications can be seen in [2], [3], 
and [13]-[24]. 

Reference [25] proposed the repetitive sampling which is 
used when there is no decision on the information of the first 
sample. This sampling found to be more efficient than the 
single sampling scheme in average sample number (ASN) 
and average run length (ARL). A rich literature is available 
on the repetitive sampling in the area of acceptance sampling 
and the control charts. The control chart based on the process 
capability index (PCI) using this sampling was considered 
by [26]. Reference [27] worked on an attribute and variable 
charts using this sampling scheme. A detailed discussion and 
applications of charts using the repetitive sampling can be 
seen in [28]-[33]. 
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According to [34] “A logic in which each proposition is 
estimated to have the percentage of truth in a subset T, the per- 
centage of indeterminacy in a subset I, and the percentage 
of falsity in a subset F, where T, I, F are defined above, 
is called Neutrosophic Logic (NL)’’. The NL is the extension 
of the traditional fuzzy approach. Therefore, the neutrosophic 
statistics (NS) is untiled the NL is the generalization of 
the classical statistics. The NS is applied when the data is 
ambiguous, vague and uncertain. Reference [35] and [36] 
introduced neutrosophic numbers in rock engineering. Ref- 
erence [37]-[40] introduced the NS in acceptance sampling 
plans. Reference [41] introduced the NS in the area of control 
chart. The variance of proposed chart using the NS was pro- 
posed by [42]. Reference [43] proposed the chart to monitor 
failure censored reliability using the NS. 

According to the best of our knowledge, there is no work 
on attribute control chart using repetitive sampling under 
the NS. In this paper, an attribute control chart using the 
repetitive sampling under the neutrosophic statistics system 
is discussed. The necessary measures of the proposed control 
chart under the neutrosophic statistics system are given. The 
control chart coefficients of the proposed control chart are 
determined using an algorithm under neutrosophic statistics 
system. The efficiency of the proposed control chart in terms 
of neutrosophic average run length (NARL) is discussed 
over the existing control chart under neutrosophic statistics 
system. From the comparison studies, it is found that the 
proposed control chart under neutrosophic statistics system is 
more sensitive in detecting a shift in the process as compared 
to the existing control chart under neutrosophic statistics sys- 
tem. An industrial application of the proposed control chart 
under neutrosophic statistics system is also given. 


Il. DESIGN OF CHART USING THE NS 

Let ny € (nj, ny} be a neutrosophic random sample is taken 
from the population having some uncertain observations, 
Dy € (Dr, Dy} be neutrosophic defective values recorded 
from ny € (nr, ny} and py € (pr. pu) be the neutrosophic 
probability of non-conforming items. Based on this infor- 
mation, [41] defined the following form of the neutrosophic 
binomial distribution 


|UCLN| n dy) ol n—dy _ . 
QN ec pO" ^N, dy 0,1... m; 


dy € (di, dy}, pw € {pL pv] | (0 
where |y| denotes the positive integer values. 

We propose the following attribute control chart using the 
repetitive sampling under the neutrosophic statistical interval 
method (NSIM). 

Step-1: Record the number of non-conforming items Dy € 
(Dr, Dy} from ny € (nr, ny}. 

Step-2: If LCLy» < Dy < UCLw»;; the process is an in- 
control state. Note here that LCLy2 € (LCLz;5, LCLy2} and 
UCLy2 € (UCL;5, UCLy2} are inner neutrosophic lower 
limits and neutrosophic upper limits, respectively. 
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Step-3: If Dy > UCLwi or Dy € LCLwi; the pro- 
cess is declared as out-of-control. Note here that LCLy; € 
(LCL; 1, LCLy,} and UCLy, € {UCL ,,, UCL y) are outer 
neutrosophic lower limits and neutrosophic upper limits, 
respectively. 

The proposed attribute control chart is the generalization 
of several control charts. For example, the proposed con- 
trol chart reduces to [41] under the NS when no repetitive 
sampling is needed. The proposed control chart also reduces 
to [27] when no uncertain observations or parameters are 
noted in the data. The neutrosophic in-control control limits 
for the proposed chart are given by 


UCLy| = npyo + kniv "Pwo (1 — puo); 
pno € {pro, puo}, kwai € tkri, kui} (2) 

UCLy2 = npyo + kn2y nPyo (1 — Pro); 
pno € {pL0, puo}, kn2 € {kr2, ku2} (3) 

npyo — Kniv"Pno (1 — Pwo); 
pno € (pro, puo}, ki € li, kui} (4) 

LCLy2 = npyo — kN2V npyo (1 — Pro); 


pno € {Pr0, puo},  kw2 € {kr2, ku2} (5) 


Note here that ky; € {kz1,kui} and ky2 € {kr2, ky2} 
are the neutrosophic coefficients. Now, we derive the nec- 
essary measures to assess the performance of the proposed 
attribute control chart when the process is in-control state 
at puo € {PLo,, Puo, }. According to the above-mentioned 
attribute control chart, the process will be out-of-control at the 
first sample when Dy > UCLy| or Dy < LCLy1. Therefore, 
the neutrosophic probability that the process is in-control, say 
pu at first sample is given by 


LCLNi 


P? yı = PILCLy| € Dy € LCLyi\pno} (6) 


Using Eq. (1), it can be written as follows 


0 |UCLNI| n dy rä 
Pim = 2 t= i 
Ni dy-Iuc p V dy Po 7 P0) 


dy —0,1,...,m dy € {dz, du}, 
PNo€ {Pr0,; Puo}, | UCLyi € {UCLy1, UCLy 4}; 
LCLyw|; € {LCL ,LCLy1} (7) 


In case, when there is no decision about the state of the pro- 
cess at the first sample, the operational process of the control 
chart will be repeated. The probability of the repetition, say 


PL is given by 


UCL n 

0 dy n—d 
Pe = J 1— N 
repN dy=|UCLy2|+1 (i ) PNo ( Pwo) 


LCLy2 n dn n—dw . 
3 250 de (5, f (1 — Pwo) ? 


UCLy| € (UCL; 1, UCLy 1}, 

LCLyw|; € {LCL11,LCLy1} : 

UCLy» € (UCL;», UCLy2}, 

LCLy» € {LCL12, LCLy2} (8) 
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TABLE 1. The NARL when ny € [50, 50]. 













































































NARLno=300 NARLno=370 
Ky [3.80810797, 0.10773763),( | [(3.4805580, 0.05542655), 
2.8508100, 0.02145566)] (3.4521577, 0.09645442)] 

Duo [0.0292,0.142] [0.0379,0.028] 

C ARLw: ARLyı 

0 [304.41, 319.69] [372.15, 374.49] 
0.01 [289.89, 290.00] [352.07, 356.48] 
0.02 [276.23, 263.46] [333.31, 339.53] 
0.05 [239.78, 199.26] [283.92, 294.37] 
0.08 [209.16, 152.59] [243.22, 256.45] 
0.1 [191.43, 128.56] [220.04, 234.53] 
0.2 (126.43, 58.58] [137.77, 154.27] 
0.3 [86.98, 29.56] [90.47, 105.73] 
0.4 [61.94, 16.22] [61.85, 75.01] 
0.5 [45.43, 9.51] [43.76, 54.82] 
0.8 [20.59, 2.53] [18.20, 24.61] 
0.9 [16.40, 1.73] [14.17, 19.56] 
0.95 [14.73, 1.44] [12.58, 17.54] 
1.0 [13.27, 1.20] (11.22, 15.78] 
1.5 [5.43, 1] [4.21, 6.40] 
2.0 [2.68, 1] [1.94, 3.15] 








The neutrosophic probability that the process is in-control 
state using the repetitive sampling is given by 


po 

0 NI. 0 0 0 

Pinn = 1 po , Pinn € [Phr Phu} (9) 
~ * reptN 


The neutrosophic average run length (NARL) when the pro- 
cess is at pyo € {Po, ; Puo,} is defined by 


1 


— T Y 
1 —€- Fin 1 
1—P reptN 


NARLyo € {NARLyz, NARLny}_ (10) 


NARLNo = 


We now derive the necessary measures when the process 
is shifted from in-control state to out-of-control state due 
to some uncountable factors. Suppose pni = (1 + )pno 
PN1 € {pxi; pui) denote the neutrosophic percent defective 
at the shifted process, where c is shift constant. Therefore, 
the neutrosophic probability that the process is in-control at 
PN1 € lpri. pui). say Po at first sample is given by 


Point = PÜLCLy: < Dy < LCLwipwo — (10) 
Using Eq. (1), it can be written as follows 


|UCLNi| n d d 
Pi = N 1 = n—ay ; 
int = È aycha (dy JP 7 Pw) 


dyn = 0, 1,..., n; dy € (dr, dy}, 
pui € [pri pui}, UCL € {UCLr1, UCLy1n)}; 
LCLn\ € {LCL11,LCLu1} (12) 


In case, when there is no decision about the state of the pro- 
cess at the first sample, the operational process of the control 
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chart will be repeated. The probability of the repetition, say 


1 . . 
P rept is given by 


UCL n 

1 dn n—d, 
P = J p 1 p N 
repN dy =|UCLn?2|+1 (h) N1 ( N1) 


LCLy2 n dy a 
SE mds (“est (1 — pi) : 


UCLy € (UCL;, UCLy1) ; 
LCLyw|; € {LCL11,LCLy1} (13) 


The neutrosophic probability that the process is in- 
control state using the repetitive sampling is given 
by 


p! 

1 Nl. 1 1 1 

Pj = a Pinn € [Phr Phu} (14) 
^ ^ reptN 

The neutrosophic average run length (NARL) when the pro- 
cess is at py € {Pii,, pui,] is defined by 


1 


p! ' 
1— za 
1-P reptN 


NARLw| € {NARLyz,, NARLuNu) (15) 


NARLyw|; = 


Suppose that roy € (ror, rou) denotes the pre-fixed values 
of NARLyo € {NARLynL, NARLny}. The values of neu- 
trosophic control chart coefficients ky; € {kz1,ky1} and 
kn2 € {kr2, ky?) and NARLn, € (NARLwr., NARLyv} are 
determined for various combinations of roy € ror, rou}, 
c and ny € {nL,ny} placed in Tables 1-3. From 
Tables 1-3, we note that values of indeterminacy inter- 
val of NARLy; €  (NARLwL,NARLywu]| decrease as 
ny € {nL,ny} increases from ny € {50,50} to ny € 
(150, 150}. For an example, when NARLyxo € 1300, 300}, 
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TABLE 2. The NARL when ny € [150, 50]. 
































































































































ARLno=300 ARLwo-370 
ky [(3.2059705, 0.23039384) | [(3.03770802, 0.03943507),( 
3.2300457, 0.07855410] 3.0061115, 0.01933103)] 
Duo [0.0162,0.0361] [0.1174, 0.2145] 

c ARLwi ARLw: 

0 [304.24, 306.77] [371.97, 388.0445] 
0.01 [286.61, 282.13] [327.69,342.437] 
0.02 [270.21, 259.80] [287.91,301.2283] 
0.05 [227.43, 204.33] [194.38,203.6801] 
0.08 [192.66, 162.37] [132.17,138.4758] 

0.1 [173.06, 140.06] [103.07,107.9297] 
0.2 [105.02, 70.96] [33.59,35.03192] 
0.3 [67.24, 39.15] [13.19,13.70168] 
0.4 [45.03, 23.17] [5.97,6.190185] 
0.5 [31.33, 14.52] [2.99,3.100434] 
0.8 [12.58, 4.61] [1,1.66] 
0.9 [9.72, 3.34] [1,1] 
0.95 [8.60, 2.86] [1,1] 
1.0 [7.64, 2.47] [1,1] 
1.5 [2.79, 0.66] [1,1] 
2.0 [1.24, 1] [1,1] 
TABLE 3. The NARL when NARLy, = 300 and 370. 
ARLwno=300 ARLwno=370 
nN [50,150] [20,200] 
kn [(3.41163588, 0.01374361),( | [(3.4033867, 0.08835755),( 
3.1956378, 0.03093748)] 3.47013388, 0.17708702)] 
pro [0.0393,0.036] [0.1324,0.0172] 

c 

0 [304.24,314.04] [370.82, 374.81] 
0.01 [287.97, 288.78] [347.78, 349.30] 
0.02 [272.76, 265.89] [326.43, 325.84] 
0.05 [232.69, 209.05] [271.18, 265.99] 
0.08 [199.73, 166.06] [226.76, 218.87] 
0.1 [180.77, 143.21] [201.97, 193.00] 
0.2 [113.72, 72.47] [117.60, 107.72] 
0.3 [75.02, 39.94] [72.38, 64.27] 
0.4 [51.51, 23.61] [46.69, 40.54] 
0.5 [36.59, 14.79] [31.34, 26.80] 
0.8 [15.38, 4.68] (11.40, 9.58] 
0.9 [12.01, 3.39] [8.53, 7.18] 
0.95 [10.68, 2.91] [7.43, 6.27] 

1.0 [9.53, 2.51] [6.50, 5.49] 
1.5 [3.61, 2.31] [2.02, 1.78] 
2.0 [1.66, 1] [1, 1] 

















ny € {50,50} and c = 0.01, the indeterminacy interval of 
NARLy| € {NARLwz, NARLyy } is [286.615, 282.1398] and 
it is [289.8948, 290.0023] when ny € (150, 150}. 

The following neutrosophic algorithm is applied to find 
the values of ky; € (Kri, kui} and kyo € {kr2, ky2} and 
NARLw| € {NARLyz, NARLwu]. 

Step-1: Pre-fix the values of ny € (ny, ny) and c. 

Step-2: Determine the values of pyo € {pxo,, Pvo.]. 
kyı € {kr1,kyı} and ky2 €  (kro,kpyo] such that 
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NARLywo € (NARLyrL, NARLNu) > row € {roL, rou}. Sev- 
eral combinations of pyo € {PL0, Puo, }, kui € {kr1, kui} 
and ky? € {kr2, ky?) will exist where NARLNo > roy. 

Step-3: Select those values of puo € lpzo.. Pvo.]. 
kny € {kr1, kui} and ky2 € {kr2, ky2} where NARLyo is 
same or very close to row € (ror, rou }- 

Step-4: Use the values of pyo € {Pxo,.Puo,}; kyı € 
{kr1,kui1} and ky2 € {kz2,ky2} to find NARLy, € 
(NARLyr, , NARLnv} for various shifts. 
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TABLE 4. The comparison in NARL when ARLy, = 370 when ny € [50, 50]. 






























































c Existing Chart Proposed Chart 

0 [871.80, 384.87] [372.15, 374.49] 
0.01 [333.48, 339.55] [352.07, 356.48] 
0.02 [298.49, 299.67] [333.31, 339.53] 
0.05 [213.04, 207.24] [283.92, 294.37] 
0.08 [152.57, 145.26] [243.22, 256.45] 
0.1 [122.81, 115.64] [220.04, 234.53] 
0.2 [45.54, 41.31] [137.77, 154.27] 
0.3 [19.77, 17.54] (90.47, 105.73] 
0.4 [9.88, 8.66] [61.85, 75.01] 
0.5 [5.59, 4.88] [43.76, 54.82] 
0.8 [1.84, 1.66] [18.20, 24.61] 
0.9 [1.49, 1.36] [14.17, 19.56] 
0.95 [1.37, 1.27] [12.58, 17.54] 
1.0 [1.28, 1.20] [11.22, 15.78] 
1.5 [1.00, 1.00] [4.21, 6.40] 
2.0 [1.00, 1.00] [1.94, 3.15] 











Ill. SIMULATION STUDY 

A control chart under the NS is said to be more efficient 
than the other if provides the smaller values of NARL at 
same parameters. The smaller the values of NARL means the 
quick indication about the shift in the process. Here, we will 
discuss the advantages of the proposed control chart over the 
neutrosophic attribute control chart proposed by [41]. 


A. COMPARISON IN THE NARL VALUES 

We first present the theoretical comparison between the pro- 
posed control chart with the attribute control chart proposed 
by [41]. For fair comparison, we set the same values of 
all specified neutrosophic parameters. The values of NARL 
when NARLyxo € {370,370} and ny € (50, 50) are given 
in Table 4. 

The Table 4 clearly indicates that the values of indetermi- 
nacy interval of NARL from the proposed control are smaller 
than the control chart proposed by [41] at all values of c. For 
an example, when c = 0.01, the values of indeterminacy 
interval in NARL are NARLy, € (333.48, 339.55} from the 
proposed control chart and indeterminacy interval in NARL 
are NARLwN; € {352.07, 356.48} from [41] control chart. 
By comparing both control charts, it can be seen that for 
the proposed control chart a shift in the process can be 
expected in between 333-339" samples while at the same 
parameters, the existing attribute control chart gives indica- 
tion between 352-356" samples. From this comparison, we 
conclude that the proposed control chart is more efficient 
than chart proposed by [41] in detecting the shift in the 
process. 


B. COMPARISON USING SIMULATION 

We also show the advantage of the proposed control chart 
over the existing neutrosophic attribute chart using the sim- 
ulated data. The data is generated from the neutrosophic 
binomial distribution with neutrosophic parameters ny € 
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FIGURE 1. The proposed chart for the simulated data. 


(150,150) and pyo € {0.1174, 0.2145}. We generated 
50 sample from the neutrosophic binomial distribution. The 
first 25 values are generated from the in-control process and 
next 25 are generated at shifted process when c = 0.3. 
For this simulation study, we have NARLywo € (370, 370}, 
ny € {150,150}, kyı € {3.0377, 0.0394} and ky2 € 
{3.0061, 0.0193}. The calculated neutrosophic control limits 
are shown in Figure 1. The values of statistic Dy € (Dr, Du} 
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FIGURE 2. The existing chart for the simulated data. 


are recorded and plotted in Figure 1. From Figure 1, it is 
clear that the proposed control chart detects shift at the 11" 
sample. The values of Dy € (Dr, Dy} are also recoded for 
the existing control and plotted on Figure 2. From Figure 2, 
we note that the existing control chart does not provide any 
signal about the shift in the process. By comparing both 
figures, we reach the conclusion that the proposed chart is 
more sensitive in the detecting the shift as compared to [41] 
control chart. 


IV. CASE STUDY 

A juice company situated in the Lahore, Pakistan is interested 
to use the proposed control chart for monitoring of non- 
conforming product. As mentioned by [44] “Frozen orange 
juice concentrate is packed in 6-oz cardboard cans. These 
cans are formed on a machine by spinning them from card- 
board stock and attaching a metal bottom panel. By inspection 
of a can, we may determine whether, when filled, it could 
possibly leak either on the side seam or around the bottom 
joint. Such a nonconforming can have an improper seal on 
either the side seam or the bottom panel” . During the inspec- 
tion, the industrial experimenter is uncertain about the classi- 
fication of some items either conforming or non-conforming. 
Due to the uncertainty, the industrial engineers can expect 
the percent non-conforming product from 0.028 to 0.0379. 
In this situation, when some observations are unclear and 
uncertain, the attribute control chart using the classical statis- 
tics cannot be applied for the monitoring of non-conforming 
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TABLE 5. Neutrosophic data of juice company [41]. 










































































Sample No.of Sample No.of 
number | nonconforming | number | nonconforming 

units D; units 

1 [12,13] 16 [8,8] 

2 [15,15] 17 [10,10] 

3 [8,10] 18 [5,8] 

4 [10,10] 19 [13,13] 

5 [44] 20 [11,13] 

6 [77] 21 [20,20] 

7 [16,16] 22 [18,20] 

8 [9.11] 23 [24,24] 

9 [14,14] 24 [15,15] 

10 [10,10] 25 [9.12] 

11 [5,8] 26 [12,12] 

12 [6,8] 27 [7.10] 

13 [17,17] 28 [13,15] 

14 [12,15] 29 [9,9] 

15 [22,22] 30 [6,9] 

items. For this real example, let ny € {50,50} and 


NARLwo € {370,370}. The neutrosophic data is taken 
from [41] and reported in Table 5. The values of Dy € 
(Dr. Dy} are plotted on Figure 3 and Figure 4 for the pro- 
posed control chart and [41] control chart, respectively. From 
Figures 3 and 4, it can be noted that some several samples are 
in the repetition and indeterminacy areas in Figure 3. We also 
note that some points are near the control limits, which needs 
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FIGURE 3. The proposed chart for the real data. 
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FIGURE 4. The existing chart for the real data. 


industrial engineers attention. Figure 4 shows that the process 
is in-control state and indicate no point near the control 
limits. 
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V. CONCLUDING REMARKS 

In this manuscript, an attribute control chart using the repet- 
itive sampling under the neutrosophic statistics system is 
discussed. The necessary measures of the proposed control 
chart under the neutrosophic statistics system are given. The 
proposed control chart is more sensitive and adequate to be 
used in uncertainty environment than the existing attribute 
control chart under classical statistics. From the simulation 
study and real example, it is concluded that the proposed 
control chart is more efficient in detecting a shift in the 
process. Therefore, the use of the proposed control chart in 
the industry will be helpful in minimize the non-conforming 
items in uncertainty. The proposed chart using other sampling 
scheme or for the big data can be considered for future 
research. 
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